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Activities under this project have included participation in the WISP program, a study of the data 
processing requirements for WISP, and theoretical studies of radio sounding, ducting, and mag- 
netoionic theory. 


WISP forticir 


Throughout the project, virtually all of the WISP team meetings were attended and the project 
PI also participated in the WISP CDR and most other aspects the planning and design efforts of 
the program. This included writing functional objectives for various WISP measurement modes 
and otherwise assisting the program PI in the planning for WISP. 


WISP Data System 


Considerable time was spent planning a data system for WISP. This system would have provided 
personal computer and data-station access to the WISP data on CD ROM and been essential in 
the analysis of the complex data that WISP would have produced. 


Sounder Theor 


An analysis of radio sounding in the magnetosphere was prepared in conjunction a NASA study 
group headed tty Dr. J. L. Green of GSFC and submitted to the NASA IMI Study Panel. One 
paper on this topic has been published, and a second manuscript on the theory and design of a 
magnetospheric radio sounder has been completed and is essentially ready for publication. This 
paper included a new theory of sounding for curved surfaces, an analysis of the effect of signal- 
to-noise ratio on sounder performance, and an analysis of the expected echoes from the magneto- 
pause, plasmapause, and plasmasphere. Copies of these papers are attached. 
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Wave Ducting 


Work on wave ducting was presented at the Kyoto General Assembly of URSI. This work in- 
cluded a comprehensive theory of ducting in different wave modes, and a manuscript on this 
topic remains in preparation. 

Magnetoionic Theory 

Although the formulation of magnetoionic theory by Stix has largely superseded that of Rateliffe 
and Budden, the Appleton-Hartree Equation is often still cited as the dispersion relation of a mag- 
netoplasma. The parameters, however, are different, and it wasn't clear how the Appleton-Hart- 
ree Equation should be expressed using the Stix parameters. The equivalent formula, in terms 
of the complex polarization ratio Ey/E,, has been found to be as follows: 


where 


PS -PL _ (PS - PL)/(P - S) - P 
P-S ~ 1 -i(E y /E x )[(P-S)ID]cosQ 


Ey = (PS -PL) sin 2 9 (PS -PL) 2 sin 4 0 + j 

E x 2PDcos0 '\ 4P 2 D 2 cos 2 0 


The usual formula results from substituting the latter into the former, and unlike the original 
Appleton-Hartree formula, this formula is applicable to all wave modes. This result new will 
be useful in future analyses of wave propagation. 
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Remote Sensing Of Substorm Dynamics 

via Radio Sounding 

P. H. Reiff, J. L. Green, R. F. Benson, D. L. Carpenter, W. Calvert, 
S. F. Fung, D. L. Gallagher, Y. Omura, B. W. Reinisch, M. F. Smith 

and W. W. L. Taylor 

ABSTRACT 

This paper describes the technique of magnetospheric radio sounding and shows how it can be 
applied to produce "images" of magnetospheric electron density distributions and their variations 
during substorms. The magnetospheric radio sounder is based on more than a half-century 
heritage of ionospheric sounding combined with modem digital techniques. Coded pulses are 
transmitted by a long dipole and the delay times and directions-of-arrival of the returning pulses are 
measured. Plasma densities from 0.1 to 10 5 cm* 3 can be remotely sensed simultaneously along 
several different directions by a digital sounder operating in the 3 kHz to 3 MHz range. Positions 
of magnetospheric plasma and magnetic boundaries, such as the plasmapause and magnetopause, 
can be monitored on a time scale of a few minutes, and plasmaspheric erosion and refilling 
quantified during a substorm cycle. Such measurements have previously been impossible to 
obtain. From a sounder suitably situated in the near-earth tail lobe, one can measure the 
magnetotail cross-section and get information on the field strength, thus allowing monitoring of 
total tail flux changes during a substorm and thus the difference between dayside and nightside 
merging rates. It may also be possible to sound a near-Earth plasmoid directly and thus sense its 
position, growth, and motion. 

INTRODUCTION 

The Earth's magnetosphere is extremely dynamic, with large-scale changes in size and shape in 
response to interplanetary conditions, and major internal reconfigurations in response to 
substorms. Single-spacecraft measurements allow determination of magnetospheric conditions 
only at specific points at any given time. Remote sensing techniques, such as auroral imaging, can 
aid our understanding tremendously by allowing monitoring of the footprint of some of the plasma 
domains. Auroral images, however, reflect only precipitating particles; particles trapped in the 
equatorial plane or those in the magnetospheric boundary layers typically do not precipitate at all. 
Remote sensing by means of far ultraviolet (FUV), extreme ultraviolet (EUV), X-ray, and 
energetic neutral atom (ENA) detection promises a great step forward in monitoring 
magnetospheric domains [Williams etaL, 1992]. 

During the NASA Space Physics Strategy-Implementation Study of 1990, we suggested using 
radio sounders to study the magnetopause and its boundary layers [Reiff, 1991], and to study the 
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plasmasphere [Green and Fung, 1994]. The magnetopause, plasm asphere, as well as the cusp and 
boundary layers, can be remotely monitored if a radio sounder is placed on a high-inclination 
polar-orbiting spacecraft with a reasonable (£6 Re) apogee (Figure 1, from Reiffetal. [1994]). 



Figure 1. Schematic contours of constant plasma density in the sunward side of the 
magnetosphere. A radio sounder on a polar-orbiting spacecraft can remotely sense the density 
structures at the magnetopause, die plasmapause, the auroral density cavity and the magnetospheric 
cusp, thus creating quantitative "images". Lower-frequency waves (dotted) probe the low-density 
surfaces, while higher-frequency waves (dashed) probe the higher-density structures. Waves on 
the boundaries yield multiple reflecting points, allowing a crude image within a few minutes and a 
composite image as illustrated within the time required for a quarter of a spacecraft orbit, by 
combining successive electron density profiles along the orbit [from Reiff et aL, 1994]. 


Radio sounding techniques date back to the study of the ionosphere by Breit and Tuve [1926]. 
Swept-frequency ground-based sounders are used to monitor the electron number density (Ne) 
structure of the ionosphere up to the F-layer density peak with high time resolution. The early 
instruments evolved into a global network which produced high-resolution ionograms (displays of 
echo delay versus frequency) on 35 mm film. The bottomside electron density profiles deduced 
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from these records provided one of the cornerstones of the success of the International 
Geophysical Year [Berkner, 1959]. 

Topside ionospheric densities, above the F-layer density maximum, were first measured in the 
Alouette/ISIS (International Satellites for Ionospheric Studies) program [Franklin and Maclean, 
1969; Jackson et al., 1980]. A series of consecutive profiles can be combined to create orbital- 
plane Ne contours from the satellite altitude to the altitude of the F layer peak density [Nelms and 
Lockwood, 1966 and Benson, 1985 and references therein]. Comparisons of densities from 
topside and bottomside sounders, and multi-spacecraft studies have indicated that the error of the 
Alouette/ISIS-derived N e values is typically within 10%, even at the most remote distances 
[Whitteker et al., 1976; Hoegy and Benson, 1988 and references therein]. ISIS demonstrated the 
synergism between sounders and optical imagers on the same spacecraft [Lui and Anger, 1973; 
Shepherd et al., 1976]. Digital technology was used in later spacecraft sounders, such as in 
Japan's ISS-b (Ionosphere Sounding Satellite), Ohzora (also called EXOS-C), and Akebono (or 
EXOS-D) spacecraft; and in the USSR's Intercosmos 19 and Cosmos 1809 missions [see Green et 
al. 1994 and references therein]. ISIS also showed the ability to remotely sound the 
magnetospheric cusp [Dyson and Winningham, 1974]. The ISIS-C spacecraft would have 
contained a sounder with a similar orbit and similar science objectives to that discussed here, but 
was never flown [J. D. Winningham, personal communication, 1994]. 

Bottomside (ground-based) ionospheric sounders have enjoyed great advances over the last 
few decades [Hunsucker, 1991]. These advanced ionospheric sounders can measure the 
frequency, time delay, amplitude, phase, Doppler shift and spread, polarization, and direction of 
arrival of the echo. Two sounders in this class are the Dynasonde (developed at the U.S. 
Department of Commerce Laboratories in Boulder) and the Digisonde (developed at the University 
of Massachusetts Lowell). These instruments offer a high degree of flexibility in measurement 
format since their operations are controlled by software [Reinisch, 1986]. Evidence that their 
scientific capabilities go far beyond the Ne profiles of the standard ionosondes is indicated from 
results that show, for example, turbulence, drifts, winds and structures [Wright and Hunsucker, 
1983; Buchau et al., 1988], A new generation of spacebome sounders can use these latest 
techniques [Reinisch et al., 1992] to make sounding of the magnetosphere possible with very 
modest power requirements by taking advantage of spread spectrum pulse compression and 
Fourier transform processing. A 500m tip-to-dp dipole is planned as transmitting antenna. 

RADIO WAVE SOUNDING OF THE MAGNETOSPHERE 

Radio wave sounding of the magnetosphere uses the same fundamental principles as 
ionospheric sounding. A cold magnetized plasma supports two freely propagating electromagnedc 
waves, the O (ordinary) and X (extraordinary) modes, each with distinct phase velocity and 
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polarization [e.g., Chen , 1974]. The propagation characteristics of these waves are determined by 
the local electron plasma frequency,)^ , and electron gyrofrequency,/$. 



Figure 2. A schematic phase velocity 
-frequency diagram showing the ordinary (O 
Mode) and extraordinary (X Mode) propa- 
gation modes supported by a cold plasma. 
(The shaded areas indicate regions of non- 
propagation). The extraordinary mode has 
two branches. The lower frequency branch 
(often called the Z mode) is trapped between 
a lower frequency cutoff f z and an upper 
resonant frequency varying between f p and 
fuh, depending on the angle of propagation 
[Ratcliffe, 1959]. The high frequency branch 
is normally referred to as the X mode. For 
the purposes of remote sensing, the O and X 
modes are used because they propagate at all 
frequencies above their cutoffs at fp and/* 
respectively. From a knowledge of the 
plasma frequency, the electron density can be 
derived; from the fx frequency, the gyro- 
frequency (and thus the field strength) can be 
derived. [Adapted from Chen , 1974.] 


Figure 2 shows the variation of the square of the ratio of the wave phase velocity v p h to the 
speed of light c on a scale of increasing frequency, for several wave modes. The cross-hatched 
regions represent the forbidden frequency ranges in which the indicated waves cannot propagate. 
The Z mode has a lower frequency cutoff f z (when v p h 2 /c 2 = *») and an upper frequency 
resonance (when Vp^/c 2 = 0) that restricts its propagation and is thus referred to as a trapped mode 
of the plasma. This wave mode is unsuitable for use in direct sounding to great distances. The 
whistler mode (not shown in the figure) is also a trapped mode. 

Since the X and O modes have no propagation restrictions above their low frequency cutoffs 
(Figure 2), they are suitable for remote radio sounding. The cutoff for the O mode is the local 
plasma frequency fp. 


r 

p 2icy e 0 m 


9 JIT, kHz 


( 1 ) 


when Ne is expressed in cm* 3 ; Eq is the permittivity of free space, e is the electron charge and m is 
the electron mass. For the X mode, the cutoff fx is 
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W4 2+ U 2+ i4 


( 2 ) 


where fg = ( l/2jt)eB/m = 28.0 B (Hz) (3) 

when B, the magnetic field strength, is expressed in nT. These two modes propagate freely with 
group velocity \ g **Vph * c when the wave frequencies are well beyond the plasma cutoffs; thus 
they are called the free-space modes. 

A swept-frequency sounder transmits and receives a sequence of X and O mode pulses, with 
the frequencies increasing stepwise. When the transmitted waves enter a region of increasing 
plasma density or magnetic field strength, the pulse is reflected when the wave frequency matches 
the cutoff frequency. Most of the inner magnetosphere has a relatively low density - thus a 
satellite lies between increasing density gradients at the magnetopause and at the plasmasphere (see 
Figure 1). For most magnetospheric applications the primary echo paths will be approximately 
perpendicular to the plasma density contours at reflection. A data record (plasmagram) consists of 
the time delay and amplitude, along with complex Doppler and angle-of-arrival information, for 
each sounding frequency and polarization. 

Thus the swept-frequency measurements permit the determination of the electron plasma 
density profiles of remote plasma regions [see, e.g., Jackson et aL, 1980; Huang and Reinisch, 
1982], so long as the density is generally increasing with distance. Magnetospheric density 
structures are quite dynamic in space and in time. Thus each transmission will typically result in 
more than one echo (Figure 3). Returns will usually be obtained from several directions for the 
same frequency transmission, but each will exhibit a different time delay, angle of return, and 
Doppler shift, thereby allowing a rough cross-section within minutes. From a sequence of profiles 
taken during a portion of a single orbit, one can produce two-dimensional, cross-sectional 
magnetospheric images (Figure 1). Raytracing studies indicate that from a dayside position much 
as in Figure 1, signal returns would also be received from the dawn and dusk magnetopause [ Fung 
et aL , 1994]. Thus information on the 3-D magnetospheric N e structures will also be obtained. 
For more information on the technique and applications of magnetospheric radio sounding, see 
Green et aL, 1993; ReijfetaL [1994], Fung et a/.[1994]. Green et a/.[1994], and Calvert et al. 
[1994]. 

THE USE OF SOUNDING IN SUBSTORM STUDIES 

A radio sounder is ideal for studying the global structure and dynamics of the plasmasphere 
and its outer boundary, the plasmapause. The plasmasphere reacts sensitively to changes in 
magnetospheric convection, with the plasmapause varying in geocentric distance between 2 and 7 
RE as magnetospheric conditions change from active to quiet [Carpenter, 1966; Spiro et al., 1981]. 
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solar wind monitoring spacecraft 
such as WIND and IMP-8, one 
can test the solar wind dependence 
of empirical magnetopause models 
such as Sibeck et al. [1991], 

Roelof and Sibeck [1993], and 
Petrinec and Russell [1993a, b\. 

In the special case of 
spacecraft near apogee in the near- 
Earth magnetotail (or a lunar-based 
instrument when the Moon is in 
the magnetotail [Reiff, 1991; 

Green and Fung, 1994]), one can 
receive the reflected signals from 
both the dawn and duskside 
equatorial magnetopause and north 
and south high latitude 
magnetopause nearly 
simultaneously (Figure 4): It can 
then determine whether the 
magnetospheric tail is flattened 
and/or twisted by the IMF, as has 
been proposed [e.g. Fairfield, 

1992 and references therein]. By 
determining the sue of the magnetotail and having some information about the field strength 
derived from the X mode echoes, it will be possible to monitor the magnetic flux in the magnetotail 
and thus the net (dayside minus tail) magnetic merging rates [Russell and McPherron, 1973]. Note 
that a sounder system located in the high-latitude lobes may be able to sound through the low- 
density plasmasheet to the far side magnetopause. 

In certain geometries, it may also be possible to monitor the growth and motion of plasmoids in 
the magnetotaiL In the near-earth neutral line model of a substorm [Hones, 1979], a plasmoid 
forms, grows, and disconnects in a substorm cycle [e.g. Baker et aL, 1987; Slavin et al, 1989]. 
If a sounder is situated in the high-latitude lobes (e.g., at Xsm ~20 Re, IZsmI -10), it may be able 
to sound off the density enhancement of the growing and moving plasmoid. A schematic of this 
reflection geometry is shown in Figure 5 (adapted from Richardson and Cowley, 1987]. In the 
near-earth magnetotail (<100 Re), the magnetotail lobe density is much smaller than the 
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Sounding tail plasmas from inside the tail 



Figure 4. A schematic of reflections from density gradients at 
the magnetopause and its boundary layers from a spacecraft 
(or a lunar base) within the cislunar magnetotail lobe. The 
concave density cavity provides an enhancement of the return 
signal and allows multiple returns from a single transmission, 
yielding a multipoint tail cross-section on less than minute 
time scales. The plasma sheet is of such low density that 
middle- (dash-dot) and high-frequency signals (dashed) can 
sound right through to reach the far boundary layer and 
magnetopause, or a low-frequency signal (dotted) can probe 
the near edge of the plasmasheet from the lobe. (Beyond 
lunar distance, most of the lobe is full of plasma and the 
density contrast at the plasma sheet is insufficient for a signal 
return, but the magnetopause should still yield a detectable 
signal). 



X -0 


-50 


-100 


plasraasheet density [Slavin et al., 1985] and a sounder 
should be able to remotely sense a plasmoid. In the distant 
(>100 Re) magnetotail, however, the lobes have a great deal 
of plasma in them and the density contrast is not sufficient 
to return an echo from the plasmoid [Frank and Paterson, 
this volume]. 

A plasmoid moving down tail past the sounder may 
allow a sequence of soundings to yield information on the 
growth and motion of the plasmoid, and monitoring of the 
deformation of the magnetopause as the plasmoid moves 
[e.g. Slavin et al., 1993]. A schematic of the possible 
geometry is shown in Figure 6. As a first test of this idea, 
we have taken the MHD plasmoid model developed by 
Omura and Green [1993] and performed raytracing 
calculations with a test sounder in the near-Earth magnetotail 
lobe. A set of sample rays emanating from the spacecraft 
are shown bouncing off of the plasmoid (Figure 7). Here 
the sounding frequency is near 1 1 kHz, corresponding to a 
density of about 1.5 cm* 3 . Divergence of the rays as they 
refract off the plasmoid indicate the defocusing aspect of the 
convex surface [Reiffet aL, 1994; Calvert et al, 1994]. 

Sounding Flasmoids 




Figure 5. A schematic of reflections c 
sounder rays (dashed) from density gradient 
associated with plasmoids. A sounde 
situated in the low-density magnetotail lob 
may be able to monitor the growth and motio 
of a plasmoid during the substorm cycl 
[adapted from Richardson and Cowley 
1987]. 



Figure 6. Schematic of sounding off a plasmoid from a cislunar magnetotail location as the 
plasmoid moves downtaiL A sounder will effectively move from locations A, to B, C, etc. The 
highest frequency transmissions (shown dashed) pass through the plasmoid and/or the plasma 
sheet and reflect from the magnetopause; middle frequencies (dash-dot) sense the plasmoid or the 
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boundary layer, whichever is closer along that line of sight The lowest frequencies (dotted) sense 
the plasma sheet or lowest-density boundary layer, whichever is nearer. 

xio* L -O Mode (at 10.55 kHz) 



x 10® 



sounder. 

Higher sounding frequencies can penetrate the plasmoid and may allow sounding of the far 
side magnetopause. Figure 8 shows a raytracing at 14 kHz (corresponding to 2.4 cm* 3 density), 
the minimum frequency which will traverse the plasmoid. Higher frequency pulses will exhibit 
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less refraction than shown here. The return signal from the far magnetopause will be weak, since 
the reflecting surface is far away; however, the signal should be focused from the concave nature 
of the magnetopause and may well be detectable. 


L-OMode (14 kHz) 



X[km] x io* 

Figure 8. Similar to figure 7, but for a sounder frequency of 14 kHz (corresponding to 2.4 
cm*3), just above the plasma frequency in the center of the plasmoid. This frequency (and 
all higher frequencies) will penetrate the plasmoid and can be used to sense the distant 
magnetopause. Higher frequency rays will be refracted less than shown here. 

CONCLUSIONS 

A magnetospheric radio sounder can provide quantitative electron density profiles 
simultaneously in several different directions. From a sequence of these profiles, contour plots of 
the density structure in the orbital plane can be constructed, with some out-of-plane information as 
well. This capability will allow the remote sensing of magnetospheric topology and plasma 
domains on minute time scales. The application of the sounder technique to substorm studies is 
straightforward and powerful. If the density contrast between the plasmoid and the background 
lobe plasma is sufficient (and it should be in cislunar space), a sounder in an appropriate location 
should be able to monitor plasmoid formation, growth and motion and thus remotely sense major 
substorm reconfigurations in the magnetotail. 
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A bstract. A radio sounder for the magnetosphere would provide nearly - 
continuous observations of the plasmapause and magnetopause, as well as cross- 
sectional density images of the plasmasphere, cusp, and auroral plasma cavity. 
Operating at frequencies of a few kilohertz to a few megahertz with 500-meter tip- 
to-tip dipole antennas and ten watts transmitter power, such a sounder would be 
capable of about 10% density resolution and 500 km spatial resolution in a few 
minutes at distances of four earth radii. A sounder for the magnetosphere is there- 
fore quite feasible and could revolutionize magnetospheric research by providing 
detailed observations of its principal density structures. 


/. INTRODUCTION 

Radio sounding involves detecting the echoes produced by wave reflections in a plasma 
at the plasma and extraordinary-mode cutoff frequencies. It thus measures density at a distance 
and is therefore ideal for measuring the density structure of the magnetosphere. 

A radio sounder uses a pulsed radio transmitter and radio receiver to measure the echo 
delay versus frequency, since this determines the density as a function of distance from the 
sounder. Operating over a suitable range of frequencies, a magnetospheric sounder could thus 
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produce continuous density measurements as a function of orbital position, yielding complete, 
two-dimensional, cross-sectional density images of the magnetopause, plasmapause, plasmasphere, 
polar cusp, and auroral cavity, as shown in Figure l 
Pioneered for the ionosphere 
during the I920's and 30's by Breit and 
Tuve [1926], Appleton [1932], Rat- 
cliffe [1933], and others, radio sound- 
ing using a world-wide network of mo- 
dem digital ionosondes has since be- 
come one of the main tools of iono- 
spheric research [Davies, 1990; Rei- 
nisch, 1986; Reinisch, et al., 1992], 

Satellite radio sounding also became 
one of the first major techniques used 
in space research during the highly- 
successful Alouette and ISIS topside- 
sounder programs of 1960's and 70's 
[Schmerling and Langille, 1969; Frank- 
lin and Maclean, 1969].. Radio sound- 
ing has thus had a long and successful Fi 8 1 Principal targets for magnetosphenc sounding 

history, and it promises equal success 

for the magnetosphere [Calvert, et al., 1972; Calvert, 1981b; Reiff, et al., 1994; Green and Fung, 
1993; Green, et al., 1993], 

A magnetospheric sounder would presumably appear as shown in Figure 2, using three 
orthogonal dipole antennas to launch radio pulses, receive their echoes, and to measure echo 
delay and direction. Antenna length is critical for radio sounding, and tip-to-tip antenna lengths 
of 500 meters for two antennas in the satellite spin plane, and 10 meters for a third antenna along 

In order to determine the feasibility of 
radio sounding it is necessary to estimate the 
echo strengths which can be produced by a 
sounder and determine the resulting density 
precision and spatial and temporal resolution. 
The characteristics of radio echoes will thus 
be discussed in the following section, fol- 
lowed by discussions of pulse transmission 
and echo reception in Sections 3 and 4, and 
an analysis precision and resolution in Section 
5. The conclusions will then be presented in 
Section 6. 


the spin axis, have been assumed for this study. 



Fig. 2 Magnetospheric sounder, using two long 
transmit/receive dipoles in the spin plane and a 
short spin-axis dipole for reception only. 
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2. ECHO CHARACTERISTICS 

The radio echoes used in sounding are produced at wave cutoffs where the refractive 
index goes to zero and total internal reflection occurs [see Ratcliffe, 1959; Budden, 1985; Stix, 
1962; Allis, et al., 1963], As a result, reflection is lossless and coherent, and the echoes are 
produced perpendicular to the refractive index contours at the point of reflection. Echo directions 
are thus determined by the density distribution of the target plasma, and are often perpendicular 
to the magnetic field at interfaces like the plasmapause and magnetopause where magnetic forces 
control the density. 


2. 1. Frequency 


A sounder launches two waves of opposite elliptical polarization known as the ordinary 
and extraordinary waves. These two waves propagate independently and produce separate echoes 
which can be distinguished by their different delay, doppler shift, and polarization. The plasma 
density can be determined from either echo, and if both are observed with sufficient precision, 
they can be used to measure the magnetic field strength at the point of reflection. 

The cutoff frequency of the ordinary mode is the plasma frequency 


f P 


1 

2 TC ^ 


e 0 m 


8,98 l N e kHz 


( 1 ) 


where N e is the electron density per Cubic centimeter, e and m are the electronic charge and mass, 
and s 0 is the permittivity of free space. The cutoff of the extraordinary mode occurs at a fre- 
quency only somewhat higher throughout most of the magnetosphere, given by 


/, - TUT* * f,n <2) 

where 

/ = = 28 B kHz (3) 

8 2itm 

is the electron cyclotron frequency and B is the magnetic field strength in microtesla. 

Since plasma densities in the magnetosphere range from less than 0. 1 cm' 3 in the auroral 
cavity and tail lobes to 10 4 cm 3 or more in the topside ionosphere, the plasma frequency in the 
magnetosphere ranges from about 3 kHz to over 1 MHz. Magnetic field strengths, on the other 
hand, range from about 30-50 microtesla in the ionosphere to about 30 nanotesla at the magneto- 
pause, and even less in the geomagnetic tail. The cyclotron frequency thus ranges from about 
15 MHz in the ionosphere to less than one kilohertz at the magnetopause, and except in the 
auroral plasma cavity [see Calvert, 1981a], it is generally less than the plasma frequency. The 
pertinent frequencies for radio sounding in the magnetosphere therefore range from a few kHz 
to a few MHz. 
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2.2 Echo Delay 

Although the phase velocities of the pertinent ordinary and extraordinary waves are greater 
than the speed of light, the group velocity v g which determines the speed of pulse propagation 
(neglecting magnetic effects) is approximately 

V, = c (1 -Ills 2 (4) 

where f is the frequency, c is the speed of light, and f p is the local plasma frequency Since the 
plasma frequency along most of the echo path is usually much less than the sounding frequency, 
the group speed is nearly the speed of light, and the round-trip echo delay is approximately 2/c 
times the distance, or about 43 milliseconds per earth radii (R E ). Echo delays in the magneto- 
sphere, at distances of up to a few earth radii, thus range from a few tenths of a second to a half 
second or more. 


2.3. Echo Power Flux 


In order to estimate echo signal strengths for a sounder, it is necessary to calculate the 
echo power flux O c (in watts per square meter) which results from a given effective radiated 
sounder power P, in the direction of the target, as follows. 


2.3.1. Planar Target. For a planar 
target, as shown in Figure 3, the power flux 
<t> x at the point of reflection is 


<!> 


P 

4 n s 2 


(5) 


where P is the effective radiated power and 
s is the target distance. Since reflection is 
coherent, this flux spreads out over a larger 
arc y = 2x when it returns to the sounder. 
Considering both dimensions of the echo 
wave front, the echo power flux <t> eo at the 
sounder is therefore 
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Fig. 3. Echo produced by a planar target where 
s is the target distance, O x is the power flux at 
reflection, and is the returning echo flux. 
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2.3.2. Curved Targe fs. In the magnetosphere, however, the target radii of curvature are 
frequently comparable to the echo distance, and it is therefore necessary to determine the echo 
power flux for curved targets, as follows. 

For a curved surface, the echo geo- 
metry becomes that in Figure 4, where s and 
s' are the object and image distances deter- 
mined by the spherical mirror equation of 
geometric optics: 


1 


_L 

s 


2 

R 


( 7 ) 


where R is the surface radius of curvature, 
positive for a convex surface like the equator- 
ial plasmapause and negative for concave sur- 
faces like the magnetopause (for a sounder in 
between these two regions). Since the arcs x 
and y in Figure 4 subtend the same angle at 
the image of the source: 



y 

X 


s + s 


=2(1 + sf R) ( 8 ) . Fig. 4. Echo for a curved target, where R is 
s the surface radius of curvature and s' is the 


Again considering both echo dimensions, the 
echo power flux thus becomes 


image distance. 


0 > = 

* 16 7ts 2 |l +5//?J|l +5/^1 


( 9 ) 


where R, and R, are the principal radii of curvature of the surface. The absolute- value bars in 
this equation are needed to keep the power flux positive for the case of a sounder outside the 
center of curvature of concave surface, since R is then negative with s greater than -R. 

Equation 9 is the equivalent of the radar equation for targets of arbitrary size. For a large 
target where R » s, it becomes Equation 6 for a planar target, whereas for a small target (R « 
s), it becomes <t> e = P(7tR,R ; )/(47ts : ) : , which is the classical radar equation for a target cross-sec- 
tion of tcR,R ; , which is the area of an ellipse with R, and R 2 as its semi-major and semi-minor 
axes. This formula applies to both regimes in the geometric approximation where the reflecting 
surfaces are smooth on the scale of a wavelength. In the large-target or "planar" approximation 
it applies to targets larger across than a Fresnel zone, given by (2 Xs)\ where s is the target 
distance and X is the wavelength. See Davies [1990] for a one-dimensional version of this for- 
mula which is applicable to ionospheric focusing. 
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Except for small-scale irregularities like those which occur in the auroral zone and at the 
plasmapause under certain conditions. Equation 9 should be a good approximation for the echo 
strengths produced by a sounder. Moreover, since such irregularities tend to produce stronger 
echoes which can also be analyzed for density [see Calvert and Schmidt, 1964], it is also a worst- 
case approximation for the minimum usable echoes produced by a sounder. 

2.4. Relative Flux 

Echo strengths for the magnetopause, plasmapause, and plasmasphere will be calculated 
by first determining <t>./P, then multiplying by the sounder power P in the direction of the target. 
O e /P will be called the relative echo flux, and it represents the effect of echo geometry, in units 
of inverse meters squared, independently of sounder power. The relative flux, which depends on- 
ly upon s, R,, and R 2 , will be calculated in this section, and P will be calculated in Section 3. 

2.4.1. Magnetopause. As shown in Figure 1, the magnetopause will be approximated by 
a 10 -R e earth-centered sphere, and the sounder radial distance will be assumed to be 6 R E . The 
echo distance and surface radii curvature are therefore s = 4 R e and R, = R 2 = -10 R E . The 
relative echo flux, according to Equation 9, is thus about 8.5 x 10 17 m'\ independent of frequen- 
cy, since the magnetopause is abrupt, and also independent of echo position, since in this approx- 
imation the sounding distance and radii of curvature are constant. 

For magnetopause distances of 8 and 12 R E , the relative fluxes become 2. 1 x 10 16 and 5.4 
x 10 17 m‘‘; which vary by less than l/s“ because of increased focusing at greater distances. Such 
focusing, which results from the s/R terms in Equation 9 for R negative and greater than s, yields 
echo fluxes about two to four times stronger than for a planar target at the same distance. As 
a result, these echoes are significantly stronger than plasmapause echoes at the same distance, 
which tend to be defocused by comparable factors. 

The actual curvature of the magnetopause will also vary with echo position, gradually de- 
creasing along the flanks of the magnetosphere. This will cause multiple echoes, and frequently 
up to five such echoes should occur, coming from front, left, right, up, and down with respect 
to the solar wind and the magnetic axis of the Earth. Multiple echoes are also expected from 
ripples on the magnetopause, sometimes producing intense focusing as the sounder passes near 
a ripple center of curvature. Such multiple echoes, which would be of great benefit in deter- 
mining the shape and behavior of the magnetopause, should be easily distinguishable from one 
another by their different delays and doppler shifts, and their directions and distances can be mea- 
sured independently, as discussed below. 

Tn addition, when the echo path misses the earth and plasmasphere, echoes should also 
sometimes be detectable from the night hemisphere, at distances of roughly 16 R s , with focusing 
gains of about 3 and relative echo fluxes of about 5 x l O' 18 m : . From a satellite at 6 R E , the 
magnetopause thus should be detectable more-or-less continuously for hours a time, and some- 
times for nearly whole orbits outside the plasmasphere. 
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2.4.2. Plasm apause. Although the ac- 
tual plasmapause may be quite different, it 
will be modeled as an L = 4 L-shell, as 
shown in Figure 5. The shape of the plas- 
mapause thus becomes 

r = L cos 2 A ( 10 ) 

where r is the radius and X is the magnetic 
latitude, and the magnetic dip angle l is given 
by 

tan/ = 2 tan X (11) 

In terms of this dip angle and the radial dis- 
tance, it can be shown that the radius of cur- 
vature in the magnetic meridian becomes 



Fig. 5. Plasmapause echo geometry, where r, 
X, and I are the radius, latitude, and dip angle; 
A and r sat are the satellite latitude and radius; 
and s and 5 are the echo distance and angle. 


*1 - r— — (12) 

(3 cos z I + l)cos / 

whereas that in the echo plane mutually per- 
pendicular to the plasmapause and magnetic meridian at the point of reflection is 


2 r cos / 
3 cos 2 / - 1 


(13) 


These two radii of curvature are shown in Figure 6 versus magnetic latitude. Notice in 
this figure that R ; becomes infinite for a magnetic latitude of 35 3°. This occurs where the dip 
angle equals the magnetic colatitude (position B in Figure 5), at the "top" of the plasmasphere, 
furthest from the magnetic equatorial plane. At this point, the azimuthal radius of curvature in- 
creases through infinity to become negative, and the plasmapause thus becomes concave in the 
azimuthal direction, like the inside of a donut hole. This causes focusing along the polar axis, 
as discussed below. The plasmapause thus changes from doubly-convex at the equator to a 
saddle-shaped surface at high latitudes, which is convex in one direction and concave in the 
other Its radii of curvature, according to Equations 12 and 13, increase from R, = 2 R E and R : 
= 4 R e at the equator to R, = 4.62 R E and R : = ± x at 35.3° magnetic latitude. Both then con- 
tinue to increase (R : increasing from minus infinity) to become R, = 5 86 R E and R : = -0 72 R E 
at 60° latitude where r = R E and the plasmapause blends into the ionosphere. 
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Fig. 6. Plasmapause radii of curvature: R, in the magnetic meridian and R 2 in the perpendic 
ular to the magnetic meridian. 


From the law of sines for the triangle at position B in Figure 5, the echo distance s, 
perpendicular to the plasmapause for a sounder at a constant geocentric altitude r sal , vanes with 
the satellite magnetic latitude A, as follows: 

» - r m (14) 

sin / 

where 5 = I + X - A is the echo deviation angle shown in Figure 5, given by 


5 = sin 1 


— sin / 

K r sat j 


(15) 


The quantities r, I, R,, and R : can thus be computed for a given magnetic latitude from Equations 
10-13, 5 then can be computed from Equation 15, and finally, s can be computed from Equation 
14. The relative echo flux can then be determined from Equation 9 and plotted against A = I + 
X - 5, as shown in Figure 7. The positions labeled A through D in Figure 5 correspond to A 
through D in this figure. 

The relative echo flux from the plasmapause for a sounder at 6 R E thus decreases from 
about 4. 1 x 10' 17 m’ ; at the equator to a minimum of 1.6 x 10 17 m \ slightly poleward of position 
B. It then increases to infinity along the axis of the plasmasphere (position C), since for a figure 
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of revolution, this is automatically where the azimuthal center of curvature always occurs. 
Although for a symmetrical plasmasphere this corresponds to echoes coming from a continuous 
ring around the pole, the predicted focusing should still occur even if the plasmasphere is not 
exactly symmetrical, at those spots near the pole where the satellite is at the local center of 
curvature. The echoes from C to D then correspond to sounding the plasmapause from the other 
side of the polar axis, at progressively lower altitudes until it merges with the ionosphere and the 
echoes disappear. The dashed curve in Figure 7 represents the same variation measured from the 
equator on the opposite hemisphere. Multiple echoes should therefore occur within about 10-20° 
of the polar axis, from high and low altitudes on opposite sides of the plasmasphere. Such 



o 30 60 90 60 30 0 


SATELLITE LATITUDE (°) 

Fig. 7 Relative echo flux from the plasmapause, where A, B, C, and D pertain to Figure 5 
and the dashed curve corresponds to a sounder in the opposite hemisphere 


2.4.3. Plasmasphere. The density contours within the plasmasphere should be different 
in shape than the plasmapause because of diffusive equilibrium. Instead of following the mag- 
netic field, according to the density model of Carpenter and Anderson [1992; see also Carpenter, 
et al., 1993], they are approximately earth-centered, with the density typically varying from about 
400 cm' 3 at the equatorial plasmapause to about 2000 cm' 3 at 2 R E geocentric altitude. The echo 
direction should thus change abruptly as the echo point moves into the plasmasphere, as shown 
in Figure l, and echoes from the plasmasphere should be detectable only within about ±60° from 
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the equator. For a satellite at 6 R E , 
the echo distance for echoes from 
the plasmasphere increases from 2 
R F to 4 R E as the radius of curvature 
decreases from 4 R E to 2 R E . As a 
consequence, for the Carpenter- 
Anderson model, the relative echo 
flux varies from about 5.4 x l O' 17 m' 

: at 180 kHz to 3.4 x 10 18 at 400 
kHz, as shown in Figure 8. 

2.4.4. Summary. For the 
three principal sounding targets in 
the magnetosphere, the estimated 
relative echo fluxes 0 5 /P, for a 
satellite at 6 R E , are as follows: For 
a nominal magnetopause at 10 R£ 
and the sounder in front of the 
Earth, the relative flux should be 
about 8.5 x 10 17 m' 2 , and should 
vary by about a factor of. 2 with 
magnetopause position. For a plasmapause extending to 4 R E at the equator, the worst-case echo 
fluxes should be about 1.6 x 10 17 m' : , and they should increase by about a factor of 2.5 from 
above the plasmasphere to the equator. Finally, for echoes from inside the plasmasphere, the 
relative fluxes should vary as shown in Figure 8. These will be used below to determine the 
nominal performance of a sounder in the magnetosphere. 

3. PULSE TRANSMISSION 

In order to determine absolute echo power fluxes, it is necessary to multiply the relative 
fluxes deduced above bv the effective radiated power of the sounder transmitter. This calculation 
will be carried out in this section by analyzing the power which can be radiated by the two spin- 
plane antennas at the pertinent frequencies 

3. 1. Equivalent Circuit 

An equivalent circuit for pulse transmission is shown in Figure 9, where is radiation 
resistance, X a is the antenna reactance, and V 3 is the root-mean-squared (rms) voltage at the root 
of the antenna. For this circuit the antenna current is 



Fig. 8. Relative flux for the plasmasphere, for the 
density model of Carpenter and Anderson. 
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Ir: + xl 


and the power delivered to the radiation resistance and hence 
radiated by the antenna is therefore 


„ . . 2 . _ A2L 

' ' ‘ “ ' Rl + Xl 


According to Jordan [1950, Ch. 14], the reactance and 
radiation resistance of a short dipole antenna, in ohms, are 
given by 


[1 - In (1/2 a)] i± 



Fig. 9. Equivalent circuit for 
transmission, where R, is the 
radiation resistance, X a and is 
the antenna reactance. 



where A, = c/f is the wavelength, L is the tip- 
to-tip antenna length, and a is its radius. 
These are shown in Figure 10, along with 
their schematic behavior at higher frequencies, 
for a 500-meter wire antenna, l mm in 
diameter Plasma effects on the antenna 


impedance and radiation for long-range 
echoes can usually be neglected, since the 
sounder operates at frequencies well above the 
local plasma and cyclotron frequencies. 


3.2. Radiated Power 


As shown in Figure 10, below about 
100 kHz the antenna reactance is quite large 
and increases toward lower frequencies as l/f, 
whereas the radiation resistance is small and 


10 6 



FREQUENCY (kHz) 

Fig. 10. Reactance and radiation resistance of 
a 500-meter wire dipole (schematic). 


decreases as f 2 . This severely limits the 

power that can be radiated at low frequencies. At 30 kHz, for example, the radiation resistance 
and reactance are approximately 0.5 Q and 9000 Q, respectively. The antenna current for V a = 
3000 volts is then only about 330 mA, and the radiated power is therefore only about 50 milli- 
watts. Although requiring careful engineering for space, high voltages are essential to radiate 
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sufficient power at low frequencies, and for such a source, according to Equations 17-19, the 
radiated power for X 3 » R. 3 becomes 

p m *X = Q-135 (20) 

1 ~ X* [1 -ln(I,/2a)] 2 

where L t is the transmitting antenna length and V max is the maximum allowed rms transmitter 
voltage at the root of the antenna. 



Fig. 1 1. Power radiated by a 10-watt, 3 kV transmitter using a 500-meter dipole antenna. 

The power radiated as a function of frequency, for L t = 500 m, a = 0.5 mm, and V max = 
3 kV, is shown in Figure 1 1, for an assumed maximum transmitter power of P max = 10 watts. 
According to Equation 20, the crossover point between voltage and power limiting occurs at the 
frequency 


ln(L/2a) - 1 *nr- (21) 

V V ^max 

which becomes about 1 12 kHz for the parameters given. This crossover frequency should be as 
low as possible for a high-performance sounder, and the 500-m antenna length was chosen to put 
this frequency down near the sounding frequencies necessary for a magnetospheric sounder. It 
varies inversely with the antenna length and the square-root of the antenna voltage, and the 
radiated power at lower frequencies varies as the fourth power of the frequency, according to 
Equation 20. 


/, = 
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Since pulse transmission is assumed to be limited by the maximum feasible root antenna 
voltage, there is no need to tune the transmitting antenna, unless this is part of the technique used 
to produce the needed high voltages. It is therefore expected that a sounder transmitter for the 
magnetosphere will consist of a high-voltage source with active circuits to reduce the voltage 
when the power limit is reached. This design also makes it simpler to cope with the impedance 
variations at the higher frequencies in Figure 1 1. 

The effective radiated power also depends on the antenna radiation pattern and on the 
polarization mismatch between the transmitted wave and the characteristic polarization of the 
wave mode producing an echo. The antenna pattern, according to Jordan [1950], yields between 
50% power gain and 25% loss within about 45° of the antenna axis normal, and the pattern null 
along the antenna axis can be eliminated by switching between the two long spin-axis antennas 
at appropriate times for the echoes being received. For the low plasma densities at the sounder, 
the pertinent wave modes are approximately circular [see Allis, et al., 1963, Ch. 4], Then, since 
the transmitted polarization of a dipole antenna is linear, polarization mismatch will reduce the 
transmitted wave signal by a factor of two. This, however, can also be compensated for by using 
both antennas at the same time. These two factors are therefore rather minor and can be 
neglected. 
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Fig. 12. Echo flux and reception thresholds for a 10-watt, 3-kV radio sounder at 6 R E , using 
500-meter transmitting and a 10-meter receiving antennas, and a 300 Hz bandwidth. 
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The expected power fluxes for the main sounding targets in the magnetosphere can then 
be calculated by multiplying the transmitted power in Figure 1 1 by the relative fluxes deduced 
above. These are shown in Figure 12, along with the reception thresholds determined in the next 
section. The scale on the right in this figure also indicates the corresponding wave field strength 
E, in volts/meter, given by E : = 3770 e , where is in W/m : . 

4 RECEPTION 


Echo reception involves receiving the pertinent echoes in the presence of external and 
internal noise, and the threshold for reception is determined by this noise. It is therefore neces- 
sary to identify the relevant noise sources as a function of frequency, as follows. 

4.1 Receiver Noise Threshold 


4. 1 I. Equivalent Circuit. The equivalent 
circuit for reception is shown in Figure 13, 
where R3 is the radiation resistance, X a is the 
antenna reactance, R, is the receiver input impe- 
dance, and L x is a series inductance which has 
been included for tuning the antenna. For this 
circuit, the received echo signal appears as a 
power P e in the radiation resistance,, given by 



where X is the wavelength, V e is the echo vol- 
tage across the radiation resistance, and X 2 /4k is 
the antenna intercept area [see Jordan, 1950, Ch. 
1 2], Without tuning (L x = 0), the echo power P r 
that reaches the receiver input terminals is then 



Fig. 13. Equivalent circuit for reception, 
where R, is the receiver input impedance 
and L x is an inductor for antenna tuning. 


p = r 

L r a 


= 


Ve 


<D e X 2 
4 % 


[ jrgq ' 

k 2 + *? 


( 23 ) 


Because of the small values of R a and large values of X a in Figure 10, tuning is essential for 
reception at low frequencies, as follows. 

4.1.2. Antenna Tuning. With antenna tuning, in which L x is adjusted to cancel X a , the 
received power, according to Equations 18 and 23, becomes 
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P. = 


* t \ 2 R a _ 5n L r $ e 


4 n R 


( 24 ) 


where L, is the receiving antenna length. Antenna tuning thus yields an echo signal power which 
is independent of frequency and increased by the large factor (X/R,) 2 . Such tuning, however, 
requires a tunable inductance of typically 50 mH for a 500-meter antenna at 30 kHz, and one of 
almost 3 henrys for a 10-meter antenna at the same frequency. Although the former should be 
relatively easy to implement with a bank of inductors and electronic switches, the latter poses an 
engineering challenge. This, however, is necessary so all three antennas can be used for 
reception to measure the echo directions. 


4 1 . 3 . Equivalent Noise Flux. The power P n of receiver thermal noise is given by 

P n = q k T ( A f (25) 

where q is the noise figure of the receiver (typically 3 or less for a good design), T is the Kelvin 
temperature of the receiver, k = 1.38 x 10' 23 J/°K is Boltzmann's constant, and Af is the receiver 
bandwidth (see Bennett [I960], or equivalent texts on electrical noise). 

It is convenient to convert this into an equivalent noise flux as if the receiver noise were 
being received by the antenna. This is accomplished by equating P„ to P r in Equation 25 and 
solving for <t> e = O n : 


5%L) 


(26) 


The thermal-noise threshold for q = 3, T, = 300° K, Af = 300 Hz, R, = 50 Q, and L r = 10 m thus 
becomes 1.2 x 10' 19 W/m 2 , independent of frequency, as shown in Figure 12. Since reception 
using the short spin-axis antenna will limit measuring echo directions out of the satellite spin 
plane, this worst-case threshold will be used to determine the signal-to-noise ratio in the 
following section It should be noted, however, that the reception thresholds will be substantially 
lower for the two longer spin-plane antennas, and hence that the spin direction should be chosen 
to favor echoes arriving approximately in the spin plane 


4 14 External Noise. The principal external noise sources are cosmic radio noise, solar 
noise bursts, terrestrial continuum, and auroral kilometric radiation (AKR). Of these, the cosmic 
radio noise is weakest and represents the minimum unavoidable noise floor for radio sounding, 
since it is steady, uniform in frequency and space, and always present. This noise, from the 
measurements of Bougeret, et al. [1984], is shown in Figure 12 for a bandwidth of 300 Hz. 
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Solar noise bursts are generated intermittently as a result of solar flares, and are often 
strong enough to preclude radio sounding when they occur Fortunately, even during solar maxi- 
mum, such bursts are sparse on time scales of an hour or more, and are therefore not considered 
a serious problem, although sounding data will be 'sst when they occur. 

The AKR, which originates in the auroral zone at frequencies of 50-700 kHz, is often 
stronger in the region where a magnetospheric sounder must operate. Also, since it is related to 
the aurora, it tends to occur at interesting times for sounding. The AKR, on the other hand, is 
invariably discrete in frequency, sometimes having bandwidths as narrow as only 5 Hz [Gumett 
and Anderson, 1981; Baumback and Calvert, 1987; Benson, et al., 1988], Unpublished studies 
of the ISEE-1 and DE-1 high-resolution wave data also suggest that there are no detectable 
emissions between its discrete components [see Calvert, 1982], It is thus feasible to eliminate 
AKR as broadcast stations are by modem digital sounders, by receiving first, then switching, if 
necessary, to an adjacent channel [Reinisch, 1986], This technique is also easy to implement, 
and over the past decade has become a standard feature of modem digital ionosondes [Reinisch, 
et al., 1992], Since this can be done quickly and there is plenty of room between the AKR 
discrete components (as shown in the figures of Gumett and Anderson), this can also be 
accomplished without significantly affecting the quality of sounding. As a consequence, AKR 
can be eliminated by frequency side-stepping and therefore should not be a significant problem 
for radio sounding in the magnetosphere 

Terrestrial continuum is produced near the equator at the plasmapause and magnetopause, 
and tends to occur relatively often at flux levels ranging from 10 :o to 10‘ 18 W/m : Hz [Kurth, 
1982; Gumett, et al., 1988; R. R. Anderson, private communication, 1994], It has also been 
found to consist of two components, as trapped continuum below the magnetopause plasma fre- 
quency in the distant tail, which is typically 30-40 kHz, and as escaping continuum above that 
frequency Escaping continuum, despite its name, is also discrete [Kurth, et al., 1981] and can 
be eliminated by frequency side-stepping, as just described. Trapped continuum, on the other 
hand, results from the time average of many slowly-drifting components and is often too smooth 
in frequency to be eliminated by this method [see Kurth, et al., 1981, Figs. 1 and 2], However, 
since by its nature it always occurs at frequencies below the maximum plasma frequencies of 
both the magnetopause and plasmapause, it should not significantly affect detecting these features 
or measuring plasma density in the plasmasphere. 

4.2. Digital Integration 

The key factor in radio sounding is detecting echoes with an adequate signal-to-noise ratio 
to measure echo directions, as discussed in Section 4.4. The signal-to-noise ratio (S/N) is the 
amplitude ratio of the received echoes to the pertinent noise threshold. The signal-to-noise ratio 
is therefore the square root of the power flux ratio between signal and noise in Figure 12. This 
signal-to-noise ratio, however, can be improved using pulse compression and spectral integration 
[Reinisch, 1986; Reinisch, et al., 1992], as follows. 


16 



Calvert, et al.. Feasibility of Magnetospheric Sounding 


(A) 


(B) 


Vv 


r~L 


MiA/i 


PHASE 




M— 

h/v 


N 


(C) 




nS 


■iVv 


^nX 


4.2. 1. Pulse Compression. Pulse compression 
involves transmitting a phase-modulated pulse 
consisting of n "chips" of length t, each chip being 
equivalent to a separate sounding pulse. The result- 
ing echoes are then phase-correlated with the original 
modulation pattern to produce a composite echo 
which is the coherent sum of the n returning echo 
chips. The accompanying noise, on the other hand, 
adds incoherently, and therefore produces a com- 
posite noise signal which is only the square root of 
n times that of a single chip. Equivalent to adding 
n repeated measurements, this increases the signal-to- 
noise ratio by the square root of n, as illustrated in 
Figure 14. This figure shows the echo (B) of a 4- 
chip phase-coded pulse (A), consisting of four 
coherent echo signals of amplitude S accompanied by 
incoherent noise of amplitude N. The four echo 
chips are then phase shifted and added (C), produc- 
ing a composite echo of amplitude nS accompanied 
by noise of amplitude n 'N. The resulting signal-to- 
noise ratio improvement is thus the square root of n 
(10 chips yielding a factor of 10’= 3.16, 100 chips 
yielding a factor of ten, etc ). 

Using suitable modulation codes which sup- 
press the contribution of adjacent echoes to the sum 
for a given echo delay, such pulse compression can 
be done for all overlapping sequences of n chips 
which have been received. It is thus possible to 
measure echoes delay in steps of a single chip, 

thereby improving the signal-to-noise ratio without sacrificing the spatial resolution. Its disad- 
vantage, however, is that it requires transmitting a long, phased-coded pulse during which echoes 
cannot be received. For example, using a 53-millisecond coded pulse consisting of sixteen 3.3- 
msec chips would yield a factor-of-four signal-to-noise ratio improvement and 3.3 msec echo- 
delay resolution at the expense of being unable to receive echoes during the first 53 milliseconds 
of echo delay. Pulse compression thus increases the dead time before echoes can be received, 
in this case precluding the reception of echoes from any closer than about (53 msec)/(43 msec/R E ) 
= 1.2 R e . Although this technique works best for echoes of long delay, like planetary radio 
echoes, it is still useful for magnetospheric sounding, using values for n of about 16 to 32, de- 
pending on the target distance. 


M 


Fig. 14. Pulse compression, in which a 
phase-coded echo (B) is correlated with 
the original chip pattern (A) to improve 
the signal-to-noise ratio by the square 
root of n. 
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Fig. 15. Spectral integration, in which m 
echo samples for the same delay are 
spectrum analyzed to improve the signal- 
to-noise ratio by the square root of m. 


4.2.2. Spectra l Integration Spectral in- 
tegration, on the other hand, involves transmitting 
multiple pulses and spectrum-analyzing the echoes 
which are received at the same delays after each 
pulse (see Figure 15 A). As with pulse com- 
pression, this yields a signal-to-noise ratio im- 
provement equal to the number of pulses which 
are summed in the spectrum analysis. The basic 
concept is that spectrum analysis using m pulses 
is equivalent to averaging m repeated measure- 
ments of coherent echoes having a specific pro- 
gressive phase shift between one pulse and the 
next. For echoes having a specific doppler shift 
or uniform phase shift caused by changing plasma 
conditions, all of the echo energy is deposited 
coherently into a single spectral channel, as shown 
in Figure 15B. Spectral integration thus produces 
m independent measurements consisting of the m 
possible echoes with different frequency offsets. 
The signal-to-noise ratio improvement for each of 
these echoes is then the square root of m, 100 
yielding 20 dB, etc. The penalty of spectral in- 
tegration, of course, is the additional time it take 
to transmit m pulses, and the trade-off thus be- 
comes temporal resolution versus spatial resolu- 
tion determined by the signal-to-noise ratio, as 
described in Section 4.4. 

4.2.3. Doppler Measurement. An extra 
advantage of spectral integration is that it also 
yields the frequency spectrum produced by the 


doppler shift of satellite motion with respect to the sounding target. For a pulse repetition period 


T p and m pulses, spectral integration yields m measurements of frequency shift over a frequency 
interval 1/T p in steps of l/mT p . For example, using T p = 0.5 sec and m = 8, it would yield eight 
echo measurements over a 2 Hz interval in steps of 0 25 Hz. Doppler shifts from zero to 2 Hz 


would thus occupy the eight channels, and negative or larger doppler shifts, if they occur, would 
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be aliased over this same 2 Hz interval. Since the velocity component v a corresponding to a give 
doppler shift Af a , is given by 


v d = 


c ^fd 

2 / 


(27) 


this would be suitable for measuring velocity components of ±1.5 km/sec at 100 kHz, without 
aliasing and with about 400 m/sec velocity resolution. For a typical satellite velocity of 1.5 
km/sec, this would permit distinguishing different echoes from eight different directions in equal 
cosine intervals from -180° to 180° with respect to the satellite velocity vector. The echo arrival 
direction can thus be determined to within about cos ' 0.25 = 75° in the forward and backward 
directions, and to within about sin' 1 0.25 = 15° at right angles to these directions. 

Although not useful for measuring echo direction, this is useful in two respects. First, it 
separates the echoes from different directions into separate spectral channels so their echo di- 
rections can be determined independently, as discussed in Section 4.4. This permits detecting 
multiple echoes from different directions at the same time, and this will be quite useful in 
distinguishing simultaneous echoes from the magnetopause and plasmasphere. Second, since the 
doppler measurement unambiguously determines the sense of the satellite velocity component in 
the echo direction, it can also be used to determine the hemisphere of echo arrival, and hence 
resolve the unavoidable 180° ambiguity of echo direction measurements from only the wave 
electric field. In other words, since it resolves this ambiguity, it permits absolute measurements 
of echo direction without magnetic antennas. 


4.2.4. Combined Pulse Compression and Spectral Integration. Using both pulse compres- 
sion and spectral integration yields an overall signal-to-noise ratio improvement equal to the 
square root of m times n: 


S/N = sfnTn SfN 0 ( 28 ) 

For example, using m = 8 and n = 16 would yield an overall signal-to-noise ratio improvement 
of V 128 = 1 1.3. This requires transmitting m phase-coded pulses of n chips each, followed by 
both lagged correlation and spectral analysis when the echoes are received. This, however, is 
easy to accomplish using modem digital computers, and the benefit, equivalent to having 11.3 2 
= 128 times more transmitter power, is worth the complexity it entails. 

4.2.5. Typical Measurement. In order to clarify how measurements will be made, the 
following specific example will illustrate how the echo data will be recorded and analyzed. An 
echo measurement at a given frequency, using 8-point spectral integration, 16-chip pulse com- 
pression, a chip length of t = 3.3 msec, and a pulse repetition period of T p = 0.5 seconds, requires 
transmitting eight 16 x 3.3 = 53-msec coded pulses in 4 seconds and sampling the echoes re- 
ceived at 3.3 msec intervals during each of the 500 - 53 = 457-msec periods between pulses. 
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This yields 457/3.3 = 135 echo samples per pulse, and 135 - 15 = 120 overlapping pulse- 
compression echo patterns for each pulse, since the last 15 echo patterns in the interval are 
incomplete. Pulse-compression thus yields 120 measurements of echo amplitude, corresponding 
to echo delays from 53 to 449 msec after each pulse ( 1 .2 to 10.4 R E ). These echo signals for the 
eight pulses are then spectrum analyzed, eight at a time for the same equal delay after each pulse. 
This requires 120 eight-point FFT calculations for each of three antennas, or a total of 360 
calculations This then produces an array 3 x 120 x 8 = 2880 complex spectral amplitudes 
representing the echoes for 3 antennas at 1 20 delays with 8 spectral shifts. The echo directions 
are then computed for each delay and spectral shift, yielding 120 x 8 = 960 measurements of 
echo amplitude and direction. The result of a sounding measurement at a given frequency is thus 
nearly a thousand measurements of echo amplitude and direction, and this process is then re- 
peated for the next frequency. The sounder, in this example, is thus capable of detecting up to 
240 echoes per second having 120 different echo delays and 8 different doppler shifts. 

Using modem microcomputers, these calculations require only a fraction of the four sec- 
onds allocated for the measurement, and can be carried out without delay while the next array 
of echo samples is being accumulated. The raw data rate, assuming eight-bit quantities are used 
for the resulting three quantities representing echo direction and amplitude, is 3 x 8 x 960 = 
23,040 bits every four seconds, or 5760 bits/sec. Using standard data-compression techniques, 
these data can be compressed by a factor of five or more, so the nominal data rate of the 
instrument is roughly 1 kbit/sec. 

4.3. Signal-tn-Noise Ratio 

The signal-to-noise ratios for the echo fluxes in Figure 12, before and after digital 
integration, are shown in Figure 16, for m = 8 and n = 16. For other integration parameters, 
these signal-to-noise ratios will increase or decrease by a factor of two for each combined factor 
of 4 for m and n; m = 4 and n = 8, for example, yielding half the indicated values, etc. 

The expected post-integration signal-to-noise ratios are a hundred or more for the main 
magnetospheric sounding targets, and this will be used below to determine the spatial resolution. 
As discussed above, the sounder necessary to achieve this requires a 10-watt, 3000-volt 
transmitter, two 500-meter transmitting antennas, and a 10-meter spin-axis receiving antenna. 
Since the signal-to-noise ratio is proportional to the square root of the echo power flux, it in- 
creases proportional to the transmitter voltage below the breakpoint frequency in Equation 21 and 
as the square root of the power above that frequency. At low frequencies, below the breakpoint, 
it increases proportional to the square of the transmitter antenna length, according to Equation 
21, and at frequencies below the noise-threshold crossover in Figure 12 it also increases 
proportional to the receiving antenna length. Although at high frequencies the only option for 
improving the signal-to-noise ratio is increasing the maximum transmitter power, at low fre- 
quencies it could be improved by increasing the root antenna voltage, antenna lengths, or both. 
These results pertain to a 300 Hz receiver bandwidth, and the signal-to-noise ratio is also 
inversely proportional to this bandwidth. 
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Fig. 1 6. Signal-to-noise ratio before (S/N 0 ) and after (S/N) sixteen-chip pulse compression 
and eight-point spectral integration. 


4.4. Direction Measurement 

Although in the ionosphere echo paths are generally vertical because of horizontal stratifi- 
cation, in the magnetosphere where echoes may arrive from any direction, it is necessary to 
measure their direction. Since these directions are approximately perpendicular to the so-called 
E-plane of the wave [see Allis, et al., 1963, Ch. 4], as shown in Figure 17, it is sufficient to 
measure these directions from three orthogonal components of the wave electric field. This re- 
quires coherent reception using three orthogonal dipole antennas, as shown in Figure 2, and this 
yields three pairs of in-phase and quadrature echo signals representing the electric-field amplitude 
and phase in three orthogonal directions. The echo direction can then be determined from the 
vector cross-product of corresponding in-phase and quadrature values, represented by E and E' 
in Figure 17, since this yields three quantities proportional to the wave amplitude and direction 
cosines of the E-plane normal. 
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Calculating echo directions thus requires a straight- 
forward vector cross product of six signals measured by 
three radio receivers, and it is best to do this on-orbit to 
conserve telemetry bandwidth. It is also stressed that the 
three antennas and phase-coherent radio receivers are not a 
luxury for a magnetospheric sounder, since all three are 
needed to measure echo directions. 

Except near a cutoff for waves almost exactly per- 
pendicular to the magnetic field, the E-plane normal 
represents an adequate approximation for the echo direction. 
Under certain circumstances, however, it will be necessary 
to correct for the wave-angle tilt x of the E-plane normal 
with respect to the wave normal, as discussed by Allis, et 
al. [1963] and shown in Figure 17. In addition, it will also 
sometimes be necessary to correct for the difference be- 
tween the wave-normal and ray-path directions, since it is 
the latter which actually determines the echo direction, and 
these can sometimes be significantly different [see Calvert, 
1966], These corrections, however, are well understood and can be determined from the local 
density, magnetic field, and wave mode determined by the sounder. The resulting ray directions 
must then be analyzed by ray-tracing and echo-distance inversion techniques which are well 
developed for the ionosphere and other applications [see Jackson, 1969a; 1969b], but remain to 
be adapted for the magnetosphere. Also, since the wave electric field is the same for opposite 
propagation directions [see Allis, et al., 1963], this leaves undetermined the sense of propagation 
along the echo direction, and it is necessary to resolve this directional ambiguity as described 
above, or from systematic changes of echo direction with orbital position. 

5. PRECISION AND RESOLUTION 

The precision of a sounder is dictated 
by the spatial, temporal, and density resolu- 
tion it can achieve, and it will be shown that 
these depend mostly on the signal-to-noise 
ratio. Other factors, like the absolute 
precision of its density and distance measure- 
ments, are determined by frequencies and 
times which can easily be measured to high 
precision with standard electronics. A 
sounder is thus automatically an absolute 
instrument, and precise calibration is not 
required. Also, by measuring the plasma 



Fig. 18 Angular uncertainty A0 of electric- 
field measurement for a signal S and random 
noise signal N. 



Fig. 17. E-plane normal deter- 
mined by the inphase and quadra- 
ture electric field vectors E and 
E', where x (usually small) is the 
wave-normal deviation angle 


22 



Calvert, et al.. Feasibility of Magnetospheric Sounding 


frequency remotely with waves, it measures the most fundamental wave property of a plasma 
without affecting the plasma being measured. As a consequence, except as discussed below, the 
only significant source of measurement error involves the ray-tracing and echo-distance inversion 
techniques mentioned above. Analyzing these errors, however, is beyond the scope of this study, 
and will be addressed separately. 

5. /. A ngular Precision and Spatial Resolution 

5.1.1. A ngular Precision. As shown in Figure 18, the electric vectors E which are used 
to determine the echo direction will be accompanied by a vector noise signal AE, equal in 
magnitude to the electric field divided by the signal-to-noise ratio. This introduces an angular 
uncertainty of ±AE/E in the electric field direction, and hence also in the echo directions deter- 
mined from E. The angular precision A0 of a sounder thus equals twice the reciprocal signal-to- 
noise ratio: 


A0=2-^ = ^- (29) 

E S/N 

A signal-to-noise ratio of 100 thus corresponding to a full-width angular precision of 0.02 
radians, or about one degree. 


5. 1.2. Spatial Resolution. As shown in 
Figure 19, the transverse spatial resolution As 
= s A0, at a distance s, is 


As = sAd = — (30) 

A S/N 

The transverse resolution for a signal-to-noise 
ratio of 100 is thus about 500 km at a dis- 
tance of four earth radii. The range resolu- 
tion, on the other hand, is determined by the 
chip length t and speed of light c: 


" 2 A/, 


( 31 ) 



where c/2 = 150 km/msec and Af t = 1/t is the Fig. 19. Transverse resolution As, determined 
bandwidth corresponding to the chip length t. by the distance s and angular resolution A0, 

A chip length of 3.3 msec thus yields 500 km and the range resolution As : determined by 
range resolution. the chip length t. 

The optimum bandwidth for receiving 

chips of length t is Af = l/t = Af„ since a wider bandwidth introduces unnecessary noise and a 
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narrower bandwidth cannot pass the necessary frequency components of the echo waveform. The 
optimum bandwidth for t = 3 3 msec is 300 Hz, as assumed for Figures 12 and 16. The sounder 
configuration introduced above is thus optimum for 1° angular resolution and 500 km spatial 
resolution in all three directions at 4 R E over the frequency range from of 30 to 300 kHz in 
Figure 16 where the signal-to-noise ratio exceeds 100. 

5 1.3. Bandwidth Selection. If the receiver bandwidth is varied, the noise power within 
this bandwidth increases proportional to the bandwidth. As a consequence, according to equation 
30, the transverse spatial resolution is proportional to the square root of the bandwidth: 

A s ± ~ yfEJ 

On the other hand, if the chip length is varied, the range resolution is proportional to chip length, 
and hence inversely proportional to the chip bandwidth: 

As„ ~ 1/A/, (33) 

Consequently, if both are varied and an optimum bandwidth (Af = Af t ) is maintained, 

As n A^ = constant ( 34 ) 

and the total volume resolution remains the same. 

This result will be useful in selecting the best bandwidth and chip length for a given 
sounder measurement. For example, if the magnetopause distance is more important than its 
transverse structure, one might use a pulse length of 3.3/5 = 0.66 msec and a bandwidth of 1500 
Hz, for a range resolution to 100 km. Equation 34 would then imply that the transverse resolu- 
tion would then increase by the square root of five, to become 1 100 km. 

5.2. Density Resolution 

The intrinsic density resolution of a sounder also depends on the bandwidth, according 
to Equation I : 

= 2 (35) 

K f 

The intrinsic resolution for a bandwidth of 300 Hz therefore varies from about 2% at 30 kHz to 
0.2% at 300 kHz. 

The actual density resolution, on the other hand, is usually determined by the frequency 
step between adjacent sounding frequencies. According to Equation 35, for example, a density 
resolution of 10% would require 5% frequency resolution, and hence l/log 10 ( 1 .05) = 47 loga- 
rithmic frequency steps per decade. 


24 



Calvert, et al., Feasibility of Magnetospheric Sounding 


5.3. Temporal Resolution 


The total time t loul for a complete measurement equals the number of frequencies times 
m times the pulse repetition period T p . For a logarithmic frequency sweep with constant frac- 
tional density resolution: 


* total m T. 


fogio (/2//1) 

p log l0 (l + A flf) 


(36) 


where f, and f : are the lower and upper frequency limits. The time required for sounding 
between 30 and 300 kHz with 10% density resolution, using m = 8 and T p = 0.5 sec, is therefore 
8 x 0 5 x I x 47 = 188 seconds or about 3 minutes. This then becomes the temporal resolution 
for measuring the magnetopause, plasmapause, and plasmasphere with 500-km spatial resolution 
in three dimensions and 10% density resolution. 


6. CONCLUSIONS 

This paper has demonstrated the feasibility of radio sounding in the magnetosphere and 
discussed the major factors affecting the design of a practical sounder... A suitable sounder would 
use three-axis dipole antennas to measure echo direction and be capable of measuring the density, 
distance, and direction of the magnetopause, plasmapause, and plasmasphere from 30 to 300 kHz, 
with 10% density resolution and 500-km spatial resolution, in about three minutes. 

The critical factors affecting sounder performance are the antenna lengths and the trans- 
mitter voltage at the root of the transmitting antennas. Using 500-meter transmitting antennas 
and a 3-kV, 10-watt transmitter, it is feasible to detect echoes at distances of four Earth radii with 
1° angular resolution. 

Using pulse compression and spectral integration, it is feasible to detect echoes at 30-300 
kHz with signal-to-noise ratios of a hundred or more, and it has been show that the angular 
precision, in radians, equals twice the reciprocal signal-to-noise ratio. It has also been shown that 
volume resolution of a sounder is constant when the chip length and bandwidth are varied. 

Multiple echoes, distinguishable by their different delays and doppler shifts, should be de- 
tected from the magnetopause, and as a result it should be possible to measure the distance and 
approximate shape of the magnetopause in a matter of seconds. The predicted strength of echoes 
from the plasmapause has also been analyzed as a function of latitude, and it has been found that 
multiple, focused echoes should occur in the polar region. A new formula has also been derived 
for specular reflection from curved surfaces, and the predicted focusing and defocusing by the 
magnetopause and plasmapause, by factors of 2 to 4, have been included in the analysis. 
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